Classical color field modified minijet model for pp and pp total cross section 
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In a recent paper, we have evaluated the gg — > gg scattering amplitude in the presence of classical 
color field generated by the colliding protons in the leading order approximation within the pQCD. 
In this work, we show that this amplitude can be resumed to obtain the classical color field modified 
gg — > gg elastic scattering amplitude. This modified amplitude is suppressed when the longitudinal 
momentum fraction, x, of the incident gluon is small. Minijet cross section is calculated using 
the modified amplitude. We show that the pp and pp cross section from \/s = 5 GeV to 30 TeV 
can be described as a sum of a hard component contributed by the modified minijet model and a 
soft component due to the exchange of the pomeron and of the 1 = exchange-degenerate to and / 
trajectories. The predicted cross section has a In 2 s asymptotic behavior which satisfies the Froissart 
bound. 

PACS numbers: 13.85.-t, 13.75.Cs, 13.60.Hb 



I. INTRODUCTION 

It is a well accepted notion that QCD is the underly- 
ing theory of hadron physics. However, computing total 
cross sections such as pp and pp cross sections in the 
large CM collision energy, y/s, remains an unresolved 
problem in QCD. Unlike hard processes which can be 
computed in the perturbation theory, to compute the 
total cross section requires knowledge of the imaginary 
part of the forward elastic scattering amplitude, which 
involves the intrinsically nonperburative zero momentum 
transfer physics. In the 60's and 70's Regge theory was 
extensively developed in attempts to understand hadron 
interactions [l|-l3jp. Even now one still cannot claim to 
have an understanding of the total cross sections from 
first principles of QCD [4j. 

We recall the minijet model was first introduced in the 
70's 0. At the time it was noticed that the rise of the 
total cross section was very similar to the jet production 
cross section. In this context, it is natural to separate 
the total cross section into a soft component and a hard 
component. Here the hard component is to be computed 
through the pQCD motivated minijet model. However, 
the minijet cross section rises too rapidly with respect to 
^/s. Since then the minijet model has been incorporated 
in the eikonal model by various authors in attempts to 
tame the rise and explain the data quantitatively [wl2|. 

In this paper, we compute of minijet cross section us- 
ing pQCD and demonstrate the taming of the rapid rise 
within the framework of QCD using the classical color 
field modified amplitude. We find the modified mini- 



jet model can qualitatively describe the total cross sec- 
tion over the entire range of the available data, i.e. from 
V« = 5 GeV to 30 TeV. 

The remainder of the paper is organized as follows. In 
Sec. UH we define the conventional minijet cross section 
and discuss its violation of the Froissart bound (FB). 
Then we present the classical color field modified minijet 
model in Sec. IIII1 In Sec. IIV1 we compare the present 
model with the data. Asymptotic behavior of the present 
model will also be discussed. This work is concluded in 
Sec. lYI Details on the derivation of the modified gg 
amplitude used in the present model can be found in [R 



II. MINIJET MODEL AND VIOLATION OF 
THE FROISSART BOUND 

Following the soft and hard decomposition of the total 
cross section we write pp and pp cross sections as (see 
also [ll|,[ll) 



Opp — (J so ft + Uhard 

<?pp= CT so / t (l + a/Vs) + (Thard 



(1) 

(2) 



For the soft component, we assume Regge theory is ap- 
plicable. Here the dominant contribution is the Pomeron 
exchange with cvp(0) = 1. So cr so ft — const x s a H )— 1 = 
const. The secondary Regge pole contribution is ex- 
pected to be dominated by I — exchange of the / and 
uj trajectories. We assume they are exchange degenerate 
trajectories [H-[l3I3 with a/(0) = 1/2 and a w (0) = 1/2. 
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1 For application of Rcggc theory in late 60's and the early 70's sec 
the two review articles [3, 0] an <3 the book Q ■ For more recent 
discussion on the Regge theory and its relation to QCD see the 
book See also those references listed therein. 



2 The idea of exchange degeneracy was first introduced in 
Soon after that, it was recognized the exchange of pairs of 
exchange-degenerate trajectories play a crucial rule in hadrons 
collision phenomenology in intermediate energy region, especially 
in the context of direct-channel and crossed-channel duality. For 
instance the presence of prominent low energy meson resonances 
in the pp channel is responsible for the presence of the 1 / y/s term 
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(Note that here the pomeron is associated with the soft 
component contribution with ap(0) = 1. This differs 
from the pomeron used in [jq . The latter is associated 
with both soft and hard contributions with an intercept 
at ap(0) = 1.08.) 

The hard component has been associated with the 
presence of jets in the final state, where the inclusive jet 
cross section is assumed to be dominated by the mini- 
jet production contribution At high energy, the 
dominant minijet cross section comes from the gg — > gg 
process and is given by: 

rl rl r—to 

<Jmnj = / . dxi I dx 2 \ di 

jHo J^L J-S+to 

Xg(xi,f4)9{x2,vl)Fmnj—jjr-, (3) 

where x\ and X2 are the longitudinal momentum frac- 
tions of the gluons from the two colliding protons. F mn j 
is introduced as a generalization factor which will be elab- 
orated in later sections. For the conventional minijet 
model, F mn j — 1. The differential cross section of the 
gluon-gluon elastic scattering at tree level is given by 

%^ = ^M>|M| 2 , where 

,,„, 9 /„ iit st us\ 9 / st us\ 

The quantities s, t and u are the Mandelstam variables 
of the subprocess that s = x\x 2 s, t — q 2 and u = — s — t, 
where q is the momentum exchange in the subprocess. 
The gluon distribution function of the incident protons 
g(x, Hg) are evaluated at the scale /i 2 = tu/s — p\, where 
Pt is the transverse momentum of the gluons in the fi- 
nal state. The parameter to is the cutoff of the squared 
exchange momentum which defines the hard scattering 
scale below which reactions are considered to be soft. 
One can verify that for large s the dominant contribu- 
tions of the gg subprocess are from the terms with t and u 
singularities. The corresponding leading order diagrams 
are the one gluon exchange in t- and tt-channcl which 
lead to the approximation in eq. ([?]). 

Violation of the Froissart bound: To examine 
the asymptotic behavior of a mn j, one can approximate 
&mnj by its most dominant contribution. Firstly, by peak 
approximate of the t and u singularities, we replace all 
the [i\ by t and ^ -> 9ira s (i ) /P . Thus, the inte- 
gral of t is asymptotically s-independent. The integrand 
of eq. (0) in the small- X region dominates since g(x, jj^) 
increases rapidly as x decreases. A conventional power 



in eq. For the pp case, the lack of low energy resonances 

leads to the absence of 1/y/s term in eq. |(TJ. For discussions on 
exchange degeneracy and its relation duality see [15H17I1 . 



law (PL) parametrization of g{x,p? g ) is g ~ x J when x 
is small. The asymptotic form of a mn j becomes 

0"mnjO) iki dx i ikt dx 2 {X1 1 2 ) J x s ' /_1 lns -( 5 ) 

This behavior is in agreement with that given in an earlier 
work (See eq. (18) of [U). The deep inelastic scattering 
data suggests J ~ 1.16 to 1.42 fl9t] so that at large s, 
a mnj violates the Froissart bound |20l [2lT ] , which requires 
that a < const, x In 2 s. 

Recently, there is a alternative parametrization for 
the gluon distribution derived directly from a Froissart 
bound satisfying fit of the proton structure function 
F2(x,Q 2 ) [22j. The distribution is a quadratic polyno- 
mial in ln(l/x) with quadratic polynomial coefficients in 
InQ 2 . Explicitly, for < x < 0.09, it reads 

xg FB (x,Q 2 ) = -0.459 - 0.143 In Q 2 - 0.0155 In 2 Q 2 

+ [0.231 + 0.00971 InQ 2 - 0.0147 In 2 Q 2 ] hx{l/x) 

+ [0.0836 + 0.06328 InQ 2 + 0.01121n 2 Q 2 ] ln 2 (l/x). (6) 

The asymptotic behavior of the minijet calculated using 
g FB (x,Q 2 )~(l/x)ln 2 (l/x)is 

fb _ f 1 , f 1 j In 2 (xi) In 2 (2:2) 6 

<J mni ~ / « x l / "^2 OC In S (7) 

JM. Jj>L Xix 2 

s as ls 

Despite the corresponding F 2 satisfying Froissart bound, 
the minijet cross section calculated with g F (x, Q 2 ) does 
not. Similar violation has been shown in neutrino- 
nucleon deep - inelastic scattering [23|. In fact, the au- 
thors in [22] pointed out that one should not expect a 
leading order approximation to be constrained by Frois- 
sart bound (HI. 

III. MINIJET IN CLASSICAL COLOR FIELD 

In the original minijet model, each incoming proton 
provides a gluon undergoing elastic scattering in vacuum 
to produce the jets. For each subprocess, the rest of the 
protons are treated as spectators. In contrast to this sim- 
ple picture, we consider the effect due to the other glu- 
ons in the protons. In particular, the large x and small 
x gluons inside the protons are treated collectively as a 
classical source and a classical color field, respectively, 
as in the color glass condensate (CGC) (for a recent re- 
view, see [25|). However, in this work, the classical field, 
instead of being a degree of freedom in the problem, is 
considered as a prescribed background of which prop- 
erties are characterized by the rapidity of the collinear 
gluons in the scattering. For a gg subprocess where the 
incident gluons have longitudinal momentum fraction x\ 
and x 2 , x\ and x 2 define the separation scales between 
the fast and slow gluons in the protons. Since x\ and x 2 
are independent variables and, to the leading order, the 
solution of the classical field is the superposition of the 
individual field generated by each proton. The effect of 
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the total classical field on the amplitude is expected to 
be a product of the individual effect. First consider the 
gluon with x\ from the proton moving along the +z direc- 
tion. The gluons in the proton with x > x\ are regarded 
as fast and treated as random classical color sources, pi, 
moving along the light-cones. The gluons with x < x\ 
are treated as a classical color field, A, generated by the 
source through the Yang-Mills equation. The field A can 
be solved in terms of pi . (The same situation appears in 
the other proton.) Therefore, the gluons are scattering 
in a background with classical color field instead of the 
vacuum. 

The classical color field interacts with the quantum 
gluons, B, involved in the hard scattering process and 
modifies the gg — > gg amplitude. Due to the color neu- 
trality assumption, the overall classical color field must 
be zero, but the fluctuation can be finite. So the lead- 
ing order contribution of the classical field to the am- 
plitude is proportional to A 2 ; therefore, B interacts at 
least twice with A. The modified amplitude depends on 
in turn, the random sources p\. To obtain a phys- 
ical amplitude, the amplitude must be averaged over the 
source with a weight function W[pi]. For that pur pose, 
we apply the Gaussian average of the IIM model [26| . Al- 
though the Gaussian average is an approximation to the 
full renormalization group solution of the weight func- 
tion, it captures the physics in both dilute and saturated 
regime. 

We refer the reader to [l3| for the detailed derivation 
of the modified gg — > gg amplitude. Here we quote the 
result of the modified amplitude. In particular, for the 
t and u-channel one gluon exchange diagrams at high 
energy, the propagator of the exchanged quantum gluon 
is modified while the vertex remained unchanged. The 
modification to the propagator due to the proton moving 
to the +z direction to the leading order in the coupling 
and the classical field contains two terms. Under the 
eikonal approximation at high energy scattering that the 
momentum exchange q^ <C v5, only one of the modifi- 
cation terms in 

.13] is relevant. The modified propagator 

is given by 

Gf u (q,x 1 )=g fl J ab ^ 

-g llI/ 5 ab ^I(Q 2 ,x 1 )sgn(q-)^ J ,. (8) 
ttL q 

The function sgn(x) — 9{x) — 6{—x) and 

r Q2 d 2 k t 



I(Q\ Xl ) 



k\ 



-p,(k 2 ± ,xi), 



where fJb{k\, x\) is given by the IIM model as 



jca s N c 



In 



1 



(9) 



(10) 



Q 2 = —q 2 and R p is the radius of proton which is set to 
be R p = 0.8768 fm. The saturation scale, Q 2 (x), char- 
acterized the density of the source and has a power law 



dependence on x. For x ~ 1, the source is dilute and 
Q 2 (x) is expected to vanish. In this analysis, we assume 

Ql{x) = Ql(^)\l-xY, (11) 

where the value of Qq, Xq and A are set to be the typical 
values from the analysis of DIS [13, HI] as Ql = 1 GeV 2 , 
xq = 3 x 10~ 4 and A = 0.28. We use p = 5 in the 
calculation since in the dilute region, Q 2 oc xg(x, Q 2 ) « 
(1 — x) p and the exponent, p — 5, is motivated by the 
simple choice of xg in 29]. In fact, the value of p is 
insensitive to the cross section since there is a very small 
contribution from the region where x ~ 1. The function 
I(Q 2 ,x) can be integrated analytically as 



I(Q 2 ,x) 



a s N c -jc 



Lio - 



Q 2 

■Li 3 (- (Q 2 s (x)R 2 p ) 



(12) 



where Li2(x) is the dilogarithm. I(Q 2 ,x) is a monotoni- 
cally increasing function in Q 2 (x) and I(Q 2 ,x = 1) = 0. 

Eq. ([5} can be further simplified. The parameter L~ 
is the average longitudinal size of the source with x for 
x% < x < 1. We defined L~ = X * P + ■ One can roughly 
estimate a reasonable value of x f° r small-x and low Q 2 



with g(x,Q 2 ). For example, using 
J = 1.2 at Q 2 = 1.8 GeV 2 , 



in Hal with 



1 



dx 



y xP^ 



f 1 x- J 

Jzi xP+ 

£ x ~ Jdx 



J -I 1 

J XiP^ 



(13) 



Therefore, x = ~ 0.17. x — 0.168 is used in 

our calculation. Moreover, when one applies the prop- 
agator to the t and u channel exchanges, the energy- 
momentum conservation of the vertex and the on-shell 
condition of the incoming and outgoing gluons requires 
q~ = q 2 1 (2xiP + ) as we show in FigJT] q~ is nega- 




FIG. 1. Momentum conservation of the vertex and on-shell 
condition for the incoming(:riP + ) and outgoing(fc) gluons: 



-2xiP -q + if 



: 



g7(2a;i P 4 



tive as in physical region q 2 < 0. Using eq. (|13p . 
q~ = q 2 L~ /(2x)- The L~ dependence in the denomi- 
nator of eq. © and in q~ are canceled. The propagator 
becomes 



1 



, a B N c 

1 2 X 7T 



I{Q 2 ,x 1 



(14) 



and has simple form of a Feynman propagator multiplied 
by a pure imaginary factor that depends on x± and Q 2 . 



4 



Iterative sum and classical field modification 
factor: The final modified propagator is given by the 
iteration of the correction in eq. (|14[) . Diagrammati- 
cally, it corresponds to the series in Fig. [5J Each blob 
represents the gluon's interaction with the A field from 
the proton and contributes the same multiplicative fac- 
tor. The sum of the series is a geometric sum. So the 
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FIG. 2. Schematic Feynman diagram which represents the 
iterative sum of the classical field modified propagator. The 
first term is the bare propagator. The blob and the two lines 
connecting to it in the series represent the interaction between 
the quantum gluon propagator and the classical field. 

final form of the propagator is the product of the bare 
propagator and a correction factor, 



tab 



.a,N r 



I -\-i^I(Q 2 , Xl ] 



2 X tt 

i(Q 2 ,xi)] + 



i-^/(QVi) 



(15) 



The expression of eq. tjl 5|) is analytic on the complex 
plane of I(Q 2 ,xi) except one point. According to an- 
alytic continuation, one can analytically continues the 
propagator to a region where I{Q 2 ,x\) is larger, there- 
fore, equivalently, the large Q 2 (xi) 2 or small X\ region. 

For the case of the other proton providing the gluon 
with X2, a similar correction factor is introduced to the 
propagator, except that the indexes — and 1 are changed 
to + and 2; namely, {L~ ,q~ ,X\] — > {L + ,q + ,x 2 }. Going 
through the same analysis for L~ and q~ , the modified 
propagator is obtained as 



G%(q, an) = g^S —f(Q , xi)f(Q 2 ,x 2 ), (16) 



where 



f{Q\x) 



l-i^I(Q 2 ,xY 



(17) 



Applying to the minijet model, the correction factor is 
multiplied to the t and u channel amplitudes as a nu- 
merical factor. It is equivalent to multiplying a factor 
F^ n Ax\,X2) to the gg — > gg differential cross section, 



F!' 



\f(Q 2 ,x 1 )f(Q 2 ,x 2 )\ 2 



(18) 



We refer to Ff£ n j as the classical color field modification 
factor. The xi-dependence of Ff~ n j is shown in Fig. [3] 



If both an and x 2 equal to one, F£ l m j — 1 and the con- 
ventional minijet is recovered. When x\ or x 2 becomes 
small, F£~ n j is l ess than 1. Therefore, the contribution 
from the high gluon density region, where x\ or x 2 is 
small, is suppressed. 




icr 5 lo- 4 o.ooi o.oi 



FIG. 3. a;i-dependence of F^ nj -,(xi, X2, Q 2 ) for X2 — 1 and 



1 GeV^ 



IV. IMPLICATION FOR pp AND pp TOTAL 
CROSS SECTION AT HIGH ENERGY 

Comparing to data: To compute the minijet cross 
section at high energy, we approximate dcr gg /dt by the 
sum of the singular terms in t and u (the last approxi- 
mation in eq. (0J). Furthermore, t and u are symmetric 
under the integral so the integral can be done by keeping 
one of them and multiplying the result by 2. The clas- 
sical color field modification is characterized by setting 
F mnj = F^ nj (xi,x 2 ,Q 2 = -i) in eq. ©. The cross 
section is calculated by 



2*0 / 2tg 



dx\ I di2 



dig(x 1 ,fj, 2 )g(x 2 ,fi 2 q ) 



-s+to 



><F,nnj(xi:X 2 ,Q 2 =-t) 

We use the running coupling 
a s (Q 2 ) 



9ira 2 (fj, 2 ) ( u§ 



4tt 



(ll-8/3)HQ 2 /A 2 QCD ) 



(19) 



(20) 



where Aqcd = 217 MeV. The hard cut-off scale is set 
to be to = 1 GeV 2 . We use the gluon distribution in 
22]. For the function I(Q 2 ,x), we use c = 4.84 as in 
the IIM model and set 7 = 0.9. So far, the minijet cross 
section includes only the gg — > gg contribution that is 
adequate for very high energy. But for smaller energy, 
one also needs to consider the contribution from qaurks 
scattering. Empirically, the ratio between the total mini- 
jet cross section and the contribution from gg — > gg can 
be approximated by 



°tot mnj 



1 



s 0.14 " 



(21) 



5 



This parameterization is motivated by the result in [30j ■ 
For the soft component of eq. (TTJ) and ([5]) , we use a so f t — 
38.5 mb and a = 1.5 GeV. 

We compare the total cross section to the pp and pp 
data for energy 5 GeV < y/s < 30 TeV. The results of 
pp and pp total cross sections from the present model 
are shown in Fig. U as the black lines with flat energy 
dependence and an initial decrease from i/s = 5 to 20 
GeV, respectively. The conventional minijet cross section 
is also shown. The figure shows that the presence of the 
classical color field suppresses the minijet cross section 
significantly. 





pp Data 
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FIG. 4. Comparison between pp and pp data with the modi- 
fied minijet model. The dash-dot line is calculated using the 
conventional minijet. The solid lines are a pp and a pv of the 
present model. The cross section of pp is higher than that of 
pp at low energy yfs = 5 to 20 GeV. The gray dashed line 
is the approximated asymptotic behavior of the total cross 
section. The data is from 31] and 321. 



One can approximate the t integral by keeping only the 
most singular term 1/t 2 . The modified minijet cross sec- 
tion is given by cr^™ = CH(s) where 



c = 



IAAitx c 



\t\ 



tHn*(\t\m 



0.489mb 



with C2(Q 2 ) being the coefficient of \n 2 (l/x) of 



g FB (x,Q 2 ) in eq. © and 



H(s) 



dx-. 



bi 2 iln 2 i 

Xl x 2 

X\X2 In 2 x\ In 2 X2 



ln 2 s 



Since our asymptotic approximation is done to the lead- 
ing order in Ins, it is correct up to a scale sq. In Fig. 
lU the asymptotic behavior of the total cross section, 
<C m (s) = Vsoft + (C/2) In 2 (s/s ) with s = 0.7 GeV 2 
is plotted to show that the numerical calculation of the 
total cross section approaches to the asymptotic form. 



Namely, a to 



> c tot cx In s. Therefore, the present 



model gives a Froissart bound saturating cross section. 
A similar asymptotic behavior were reported recently in 

fum. 



Suppressed in the small x region: At the tree 
level, the minijet is motivated by the exchange of gluon 
with Q 2 . This exchange gluon can be considered as 
a probe of the protons. When the scattering involves 
collinear gluons at small xs, many partons with a higher 
x value become a part of the source and the correspond- 
ing classical field increases. The exchange gluon has to 
interact with the strong classical field before reaching the 
collinear gluon. Throughout the process, the propagator 
acquires a suppressive factor which depends on x. 



V. CONCLUSIONS 



Note that, instead of optimizing the quality of the fit, 
our calculation is intended to show that, with the con- 
sideration of the effect of the classical color field from 
the high energy protons, the rise of the cross section can 
be satisfactorily described by the present model with a 
reasonable choice of parameters. 

Asymptotic behavior of the cr mn j-. Beside the 
ability to fit the data, the present model exhibits a Frois- 
sart bound satisfying asymptotic behavior. Consider the 
integral of eq. (|I9[) in the small x\ and xi region where 
Q 2 (x) becomes large. The dilogarithm Li2(—z) in eq. 
(IT2l) is proportional to ln 2 (z) for large \z\. So that the 
power law behavior of Q 2 (x) leads to a logarithmic be- 
havior of I(Q 2 ,x) cx ln(x) for small x. Thus, Ffil n j, for 
small x, is 



A 2 In 2 (~iR 2 ) In 2 xi X 2 In 2 {-tR 2 ) In 2 x 2 



The minijet model was originally introduced based on 
a striking resemblance between the rise of the total cross 
section with energy and the rise of jet (dominated by 
minijet) cross section. In this work, we show that within 
the framework of QCD, the interaction between the gluon 
and the classical color field can play an important role in 
understanding this rise. 

Among different paramctcrizations of the gluon distri- 
bution function, we choose the one that is consistent to 
the Froissart bound satisfying F 2 derived in [22| . We first 
show that, even with this distribution, the conventional 
minijet cross section still violates the Froissart bound. 

We then apply the formalism of the gluon-classical 
color field interaction within the Gaussian approximation 
of CGC. In this formalism, to leading order in the cou- 
pling constant and the strength of the sources, we found 
that the classical color field introduces a modification 
factor F£° n j in the integrand of the minijet cross section 
that leads to a suppression of the rate of the rise of the 
total cross section. This allows a good description of the 
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rise of the cross section comparing to the data as show in 
Figure 2J In addition, the total cross section has a In 2 s 
asymptotic behavior which satisfies the Froissart bound. 
In this picture, the rise of the minijet cross section is still 
driven by the growth of gluon density in small x, but it 
is suppressed by the quantum-classical interaction from 
the dense medium. 

Although the complete solution of the classical color 
field due to two colliding protons has not yet found, our 
approximated solution to the leading order in the source, 
p and coupling g with Gaussian average illustrates a sig- 



nificant effect due to the quantum-classical interactions. 

In the present model, the collinear factorization is as- 
sumed. In general, including the interaction with the 
classical field from the approaching protons would break 
the collinear factorization. However, in [13j . we found 
that the Feynman diagrams consisting of interactions be- 
tween the collinear gluon and the classical field vanish. 
The only correction to the gg gg amplitude is of the 
form of eq. (jTHJ) which is factorisable. This justifies the 
factorization assumption. Nevertheless, this factoriza- 
tion property is not guaranteed in the higher order cal- 
culation. 
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